o 

(N 



6 



r- 

^ 



S 



Embryo impacts and gas giant mergers I: Dichotomy of Jupiter 

and Saturn's core mass 



2,4 



2 



Shu Lin Li^'^, C.B. Agnor^, and D. N. C. Lin 

^Department of Astronomy, '^Kavli Institute of Astronomy and Astrophysics, Peking 

University, Beijing, China 

Department of Astronomy and Astrophysics, University of California Santa Cruz, USA 



, ^Astronomy Unit, School of Mathematical Sciences, Queen Mary University of London, 

^ ! United Kingdom 



ABSTRACT 

Interior to the gaseous envelopes of Saturn, Uranus, and Neptune, there are 
high-density cores with masses larger than 10 Earth masses. According to the 
^ . conventional sequential accretion hypothesis, such massive cores are needed for 

the onset of efficient accretion of their gaseous envelopes. However, Jupiter's 
gaseous envelope is more massive and core may be less massive than those of 
^ ! Saturn. In order to account for this structural diversity and the super-solar 

CN ! metallicity in the envelope of Jupiter and Saturn, we investigate the possibility 



that they may have either merged with other gas giants or consumed several 



r-^ , Earth-mass proto-planetary embryos during or after the rapid accretion of their 

O \ envelope. In general, impinging sub-Earth-mass planetesimals disintegrate in 

gas giants' envelopes deposit heavy elements well outside the cores and locally 
suppress the convection. Consequently, their fragments sediment to promote the 
k> \ growth of cores. Through a series of numerical simulations, we show that it is 

possible for colliding super-Earth-mass embryos to reach the cores of gas giants. 
Direct parabolic collisions also lead to the coalescence of gas giants and merging of 
their cores. In these cases, the energy released from the impact leads to vigorous 
convective motion throughout the envelope and the erosion of the cores. This 
dichotomy contributes to the observed dispersion in the internal structure and 
atmospheric composition between Jupiter and Saturn and other gas giant planets 
and elsewhere. 

Subject headings: planetary systems 
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Introduction 



I n the conventional se quential accretion hypothesis (SAH) for gas giant planet formation (jSafronov 



19691 : iPoUack et al.l Il996[ ). heavy elements first coagulate into planetesimals and proto- 



planetary embryos. When they attain a critical mass MpHt, ^ 3 — lOMm, they begin to 



accrete gas efficiently and to evolve into gas giants (IStevensonlll982l : iBodenheimer &: Pollack 



19861). 



Perhaps the strongest supports for SAH are the inferred existence of so hd cores and 
super- solar metallicity in the envelopes of solar-system gas and ice giants ( jGuillot et al. 



2004J : iMilitzer et al.l 120081 ) . These models are constructed to match the observed masses, 



radii, and moments of inertial of present-day Jupiter, Saturn, Uranus, and Neptune. 

Despite this verification, theoretical explanations are needed for two observed properties: 
1) a significant difference in the core mass between Jupiter and Saturn and 2) the super- 
solar metallicity in the envelopes of gas giants. In §2, we discuss the uncertainties and pose 
challenges to the most simple version of SAH introduced by these two puzzles. We suggest 
that the dispersion in the core mass and the enhanced metallicity in the envelope of gas 
giants may be the consequence of some catastrophic coUisional events which occurred well 
after the onset of efficient gas accretion. We discuss the common occurrence of giant impact 
and merger (GIM) events during the formation of isolated and multiple proto gas giants as 
well as throughout their long-term dynamical evolution. 

We carry out numerical simulation to demonstrate that GIM's can indeed lead to diverse 
planetary structures. In §3, we briefly describe our numerical methods and discuss the 
typical range of appropriate orbital parameters and planetary properties. With a smoothed 
particle hydrodynamic (SPH) scheme, we illustrate the possibility that sufficiently massive 
protoplanetary embryos may be able to intrude through the massive envelope and impact 
onto the cores of some gas giants. With three examples of head-on, parabolic collisions 
between a Saturn-like gas giant with super-Earth embryos or another gas giant, we show 
that spherical symmetry is quickly restored after the the cores of the merger products are 
thoroughly mixed. 

In order to extend the simulation to a larger dynamical range in density and over 
protracted phases of dynamical evolution, we construct an spherically symmetric Lagrangian 
hydrodynamic (LHD) scheme. With this LHD scheme, we consider, in §4, impactors with 
a range of masses and head-on collisional speeds. We show that low- velocity and low-mass 
impactors generally disintegrated in the gaseous envelope of gas giants and their debris has 
a tendency to sediment towards the planets' cores. Such events lead to the enlargement 
of the core mass. However, large-mass and high-speed impactors can survive their passage 
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through the gaseous envelope and reach the gas giants' core. The dissipation of kinetic 
energy into internal energy induces efficient and vigorous thermal convection and erosion of 
the cores. These diverse outcomes may be responsible for i ntroducing a v ariety of internal 
structure for long-period gas giants. In a follow-up paper II (JLi et al.ll2010[ ). we will apply a 
similar scenario to a series of models for short period planets with diverse internal structures. 
Finally, we summarize our results in §5. 



2. The diverse structure of gas giants due to giant impacts and mergers 

during proto-planetary formation 

In the sequential accretion hypothesis collisions represent a principal mode of growth 
and planetary evolution. Giant planetary collisions have been suggested to explain a variety 
of bulk planetary characteristics in the solar system , including the origin of Earth's Moon 
(JHartmann fc Davis! Il975t ICameron fc Ward! Il976l ) , collisional stripp ing of Mercury's pr i- 



mordial mantle ( jBenz et al. 



1988 



20071 ) and Uranus' large obliquity ( ISlattery et al.lll992l ). 



Giant impacts during the accretion epoch display a diversity of outcome morpholo gies, each 



with different implications for the development of planetary characteristics (e.g. IWetherill 
19851 : lAgnor &: Asphaugl l2004l : lAsphaug et al.ll2006l ). Recent studies have begun to explore 
how these energeti c impacts influence the luminosity, detectability and characteristics of 
extrasolar planets ( lAnic et al.l 120071 : llkoma et al.ll2006l ). Here we build on these works and 
examine how giant planetary collisions may account for the diversity of the interior struc- 
ture of gas giant planets, and the differences in core mass between Jupiter and Saturn in 
particular. 

In this section, we provide arguments to highlight the critical roles of giant impacts and 
mergers during the formation and evolution of gas giant planets and the potential outcomes 
of such events. 



2.1. Observed structural diversity and its paradoxical implications 

Although they are both classified as gas giants, Jupiter and Saturn have very different 
internal structures. Within considerable uncertainty, the estimated core mass of Saturn is 
more massive (by up to three times) than that of Jupiter. In contrast, Jupiter's total mass 
(which is mostly contained in its envelope) is around three times that of Saturn. The average 
metallicity of both planets are inferred to be several times that of the Sun, albeit they are 
unlikely to be identical. 
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According to SAH, gas giants acquired their cores prior to the accretion of their gaseous 
envelopes. Current quantitative models indicate that a critical co re mass (Mrrit^ lO Mm) is 
needed for the onset of efficient accretion of the surr ounding gas (IPqllack et al.lll996l ). This 



mass is determined by the heat transfer effic i encv (llnaba fc Ikoma 



boundary condition s ( jBodenheimer et al.ll2000[ IWuchterl et al. 
20061 : iRafikovlbooeh . 



I2OOOI: 



20031) and the nebula 



Papaloizou fc Terquem 



There are considerable uncertainties on the magnitude of Merit- The first sets of quan- 
titative models were constructed under the assumptions of spherical symmetry, hydrostatic 
equilibrium, ef ficient convection, in terstellar-grain opacity, and idealized planetesimal bom- 
bardment rate (jPoUack et al.lll996l ). These models indicate that uninterrupted planetesimal 
bombardments and a cooling barrier can delay the formation of Jupiter in a minimum mass 
nebula (MMN) by 10-20 Myr. 

However, gas giants can only acquire their massive envelope inside gas-rich protostellar 
disks, i.e. their formation time scale Tm must be smaller than the typical gaseous disk 
depletion time scale rdep(~ 3 — 5Myr). There have been several previous attempts to resolve 
this theoretical challenge. Alibert et al. (2004) suggested that type I migration of Earth- 
mass embryos (Ward 1997) may speed up the time scale needed for the onset of efficient 
gas accretion. Although this scenario may provide a solution for Saturn (because its core is 
relatively massive), it may not be appropriate for Jupiter if its core mass is limited to a few 



In an attempt to account for Jupiter's low- mass core, llnaba fc Ikomal (120031 ) : lHubickyj et al 



(120051 ) constructed models to show that the gas accretion rate increases with the heat trans- 
port efficiency. Since there are radiative regions in proto gas giants' envelope, the critical 
mass for the onset of efficient gas accretion increases with the opacity of the infalling enve- 
lope. The main sources of opacity in these regions are grains. Grain growth and formation 
of planetesimals can reduce the opacity of residual disk gas substantially (see §2.2p . It is also 
possible for dust to coagulate and sediment from the tenuous outer envelope of proto-planets 
jHelled et allboosh . 



In principle, these effects can suppress the barrier for the onset of gas accretion. How- 
ever, a rich populati on of super-Earths (with masses Mp > a few M m) has been discovered 
with radial- velocity ( JMayor et al.ll2009l ) and transit ( JQueloz et al.l 12009) surveys. Some of 
these super-Earths have masses larger than Jupiter's estimated core mass (at least for some 
models). This discrepancy would introduce a paradox, if we assume these super Earth at- 
tained their present-day mass and became "failed cores" in gas-rich protostellar disks because 
their mass Mp < Mcnt- It is possible that the onset of gas accretion onto these super-Earths 
may be delayed by planetesimal bombardment, albeit such events are likely to be suppressed 
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by the formation of a gap in the planetesimal disk (see §2.41 IZhou fc LinI 120071 ) 



Although this accretion time scale issue is not the central focus of our discussion here, 
scenarios proposed for its resolution have implications which exacerbate the second paradox, 
i. e. gas giants' structural diversity. If the internal structure of the mature planets is assumed 
to be the result of some critical conditions that led to their formation, the large observed 
dispersion would imply the existence of many possible paths for their final assemblage. 

For example, both Jupiter and Saturn probably formed in a similar solar nebula en- 
vironment close to their present-day location (albeit at different times). It is natural to 
assume a similar magnitude for Merit for the onset of efficient gas accretion, especially 
i f their separation may have widen through planetesimal scattering since their formation 
( iHahn &: Malhotra] Il999l ). The dichotomy of their internal structure is indicative that it 
only partially reflects the critical initial condition for the onset of efficient gas accretion. It 
is entirely possible that the present-day internal structure of gas giants has been affected 
by some local or temporal stochastic processes after they have initiated (or terminated) the 
process of rapid accretion of their massive envelope. 



2.2. Heavy elemental abundance in gas giants' envelope 

A related conundrum is the super-solar metallicity content of Jupiter and Saturn's 
gaseous envelope. In the SAH, heavy elements are assumed to be mostly condensed into 
grains, efficiently segregated from the disk gas, and converted into planetesimals, embryos, 
and cores prior to the onset of the gas accretion. The abundance distribution of refractory 
elements in the meteorites suggests that the solar nebula had an initial composition similar 
to that of the sun. The emergence of solid cores with Mp > Merit would decrease the nebula's 
overall metallicity, including contribution from molecules and grains. If gas giants' envelopes 
are composed solely of nebula gas and grains accreted by their cores, they would attain a sub- 
solar metallicity. This inference clearly contradicts with the observed super-solar metallicity 
in Jupiter and Saturn's envelopes. 

One possible solution for this paradox is that Jupiter and Saturn formed in regions 
where heavy elements are locally enhanced. Due to hydrodynamic drag, grains formed in 
the outer disk regions migrate inward. They may stall and accumulate near the snow line 
(flico where the disk mid-plane temperature decreases belo w ice condensation temperature at 
160K, see e.g. ICuzzi fc Zahnld 12004 : iKretke fc Linll2007| ). This migr ation barrier to g rains 
and embryos provides a preferred location for gas giants to assemble (jlda &: Linll2008bl ) (see 
§2.3p . However, a strong metallicity enhancement would increase the effective opacity of 
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gaseous envelopes around solid cores and the magni tude of Mnri+ jinaba fc Ikomal 120031) we ll 



beyond the current estimated core mass for Jupiter (JGuillot et al.ll2004l : iMilitzer et al 



20081 ) ■ 



Although coagulation and sedimentation of grains in protoplan ets' envelopes may provide 
a potential solution for Jupiter's small core ( iHelled et al.ll2008[ ). it is not clear why such a 
process did not also enable Saturn to attain a small Merit and low core mass. 

Another possible cause of the high metallicity in gas giants' envelopes is the erosion 
of their cores. However, such a scenario would require much more massive cores (up to 
~ 30 — SOM^) prior to the onset of efficient gas accretion. In a MMN, formation of such 
massive embryos is challenging though not impossible. Core erosion by convective gas may 
also be suppressed by the compositional gradient in the boundary layer separating gas giants' 
cores and envelopes (JGuillot et al.ll2004l ). This scenario can also be easily ruled out by any 
dip in the planetary mass distribution in the range of 10 — lOOM^ (jida &: Lirul2004l ). 



2.3. Emergence of cores of the first-generation gas giants 

In this paper, we propose giant impacts of residual embryos onto gaseous planets during 
and after the advanced stages of their envelope accretion as a potential mechanism which 
may lead to the resolution of the structural diversity (see §2.ip and metal-rich envelope (see 
§2.2p . In order to distinguish between the ab initial core and the giant-impact scenarios, 
we first discuss the dynamical evolution of massive embryos in the context of gas giant 
formation. 



The growth of embryos is hr nited by their iso lation mass ((Mjso, see e.g. iLissauerl 119871 ) 
and long growth timescale (xm) (jIda fc Lirul2004l ). In a MMN- type protostellar disks, the 
most massive embryos emerge beyond Oice within the gas depletion timescale (rdep) is ~ 3 — 
5M® < Merit. Furthermore, the onset of gas giant formation is challenged by the type I orbital 
migration of their super-Earth embryonic cores. The rate and direction of type I migration 
is determined by both the surfac e density (S„) and temperature T (rn ore accurately the 
entropy) gradient in the disk gas (JMasset et al.ll2006l : |Paardekooperll2010[ ). These quantities 
are determined by the efficiency of turbulent-induced angular momentum transport. Po ssible 
causes of t urbulence include the magneto- rotational (MR I) (JBalbus fc Hawleyl Il99ll ) and 
convective (JLin fc Papaloizoulll980l : iLesur fc Ogilviell2010l ) instabilities. 



Near the snow line of a MMN, the ionization fraction near the disk midplane is negligible 
and the gas there is unaffected by MRI turbulence. However, the ionization fraction may 
be much higher near the disk su rface layers wh ich are exposed to the ionizing photons from 
the central stars and elsewhere (jGammielll996[ ). The depth of layers which are affected by 
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the MRI turbulence is determined by an ionization equilibrium in which the ionization rate 
is balanced by the rate of recombination on the grains. Interior to a^ce, the sublimation 
of volatile ices eliminates a majority of the recombination sites, enhances the ionization 
fraction of the disk gas, and increases the depth of the surface layer which is regulated by 
MRI turbulence. Across aice, effective viscosity has a negative gradient, which in quasi steady 
state, leads to a positive gradient in Eg and Pg. For a range of disk mass accretion rates 
(Mrf ~ 10~^ — 10~*Mq yr~^), the local Pg maximur n at Oinp provides a n atural barrier for the 
orbital decay of grains against hydrodynamic drag (JKretke fc Liru 120071 ). This structure also 



rever ses the direction of type I migration due to the effect of horseshoe drag (JMasset et al. 



20061 ). O utward migration o f embryos by a few AU may also be induced by an entropy 
gradient (|Paardekooperll2010l ). 



Thus, similar to grains' orbital evolution, embryos' migration leads to them to congre- 
gate near Oice where both the gas surface density (S^) and pressure (Pg) attain local maxima 
in the disk due to the transition in the efficient turbulent angular momentum transport. 
Through the accumulation of neighboring embryos and planetesimals a protoplanet may 
acquire sufficient mass {i.e. ~ Merit) to evolve into a core which efficiently accretes gas 
Jida fc LinlboOSbh . 



2.4. Embryos' dynamical evolution around isolated gas giants 

During the early phase of gas accretion, the time scale for a proto-planet to double its 
mass is longer than the synodic periods (relative to its orbit) of planetesimals and embryos 
within its feeding zone. Its tidal perturbation on the nearby embryos and their interaction 
with the disk gas induces their orbit to e volve away from that of the dominant proto-planet 
(IZhou fc LinI 120071 : IShiraishi &: Idall2008[ ). Therefore slowly ( "adiabatically" ) growing proto- 
planets cannot accrete those embryos within their expanding feeding zones in contrast to the 
assumption made in previous models (jPoUack et al.lll996l : iDodson- Robinson et al.l 120081 ). 



This initial suppression of giant impacts reduces the energy released by the growing 
proto-planet and promotes the gas accretion rate. In the limit of negligibly small thermal 
feedback and weak proto-planetary torque, proto-planets' uninhibited gas accretion rate 
(Mp) may be approximated by the Bondi formula such that the mass doubling time scale 
would be 
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where H is the disk thickness. 



The above expression indicates that, once it is initiated, Bondi-type gas accretion is 



a runaway process. During their growth, there are several effects which can lead to GIM 
events: 

Al) Proto gas giants' accretion rate Mp is eventually suppressed by the structural changes 
in their natal disks induced by their tidal perturb ation when their M„ excee ds ~ {H/aYM^, 
which is comparable to Mj in a MMN near aice (JDobbs-Dixon et al.l 120071 ) . However, just 
before they have acquired an a dequate mass to induce gap formation, proto-planets' r^ 
reduces below a thousand years (JTanigawa fc Watanabell2002l ). which is comparable to the 
synodic periods of their nearby planetesimals and embryos (with semi major axis a ± 6a ^ 
0.1 — 0.2a). Gas giants' rapid (non adiabatic) mass increases destabilize these embryos' orbits 
and a modest fraction of these objects collide with and are captured by the emerging gas 
giants Jzhou fc LinlboOTh . 



A2) After gas giants' accretion is quenched by gap formation, their gravity perturbs the 
orbital elements of residual planetesimals and embryos out to quite large distances. To 
the lowest order, gas giants' secular perturbation leads to the modulation of the embryos' 
eccentricity and the precession of their longitude of periastron. The am plitude of embryos' 
eccen tricity modulation increases with gas giants' eccentricity Cp and Mp (JMurray fc Holman 
19991 ). Many known gas giants have significant e^'s and they can induce modest amplitude 
eccentricity modulations among their nearby embryos. 

Embryos' orbital precession frequency is determined by their a's. Due to the self gravity 
of their natal disks, orbits of gas giants and residual embryos also precess with frequencies 
determined by their a's relative to the T,g + T,j_ distribution of the disk. Secular resonances 
occur at special locations where the embryos' precession frequencies due to the gravitational 
perturbation of the gas giants matches their differential precession frequencies due to the disk 
potential. In these resonances, embryos' eccentricities can be greatly excited by the nearly 
constant differential longitude of periastrons. During the depletion of the disk, gas giants' 
secular resonances sweep over w ide regions a nd destabilize the orb i ts of the planetesimal s 
and embryos at these locations (IWardlll98ll : iNagasawa et al.l l2005l : iThommes et al.l l2008l ) . 
Some embryos may cross the orbits of and collide with the gas giants. 



2.5. GIMs in systems with multiple gas giants. 

There are two gas giants, Jupiter and Saturn, in the solar system. Several stars (such as 
V And) have known multiple planets around them. Around other planet-be aring stars, there 



are in dications of additional planetary companions with comparable M„'s (jCumming et al. 



20081 ). In these systems with multiple planets, there are several additional avenues which 
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may lead to catastrophic impacts of residual embryos onto their envelopes. For example: 

A3) Regardless of the low formation and survival probability of the first-generation planets 
( jlda fc LinI l2008bl ) . their birth leads to gap fo rmation which provide a secondary pressure 



maxi mum in the gas disk outside their orbits (JLin fc Papaloizoul Il979l Il993l : iBryden et al. 



2000l ). This induced disk structure also sets up a barrier against the orbital decay of grains 
and embryos. The accumulation of heavy elements beyond the outer edge of the gap promotes 
the formation of more-distant second-gener ation cores ori time scales much shorter than that 
required for them to assemble in isolation ( jBryden et al.ll2000l ). In this sequential formation 
scenario, the second-generation gas giants' rapid mass increase, during their advanced growth 
phase, also introduces non-adiabatic changes in the systems' gravitational potential which 
have a tendency to destabilize the orbits of the nearby residual embryos. 

A4) After the gas depletion, planets' orbits continue to evolve due to their scattering of 
residual planetesimals. If multiple planets' orbits have spread out from a more compact 
configuration, orbital instabilities and sc attering amongst gian t planets may also lead to late 
accretion of planetesimals by gas giants (jTsiganis et al.ll2005[ ). 



2.6. Disruption of impacting embryos 

The main focus of this paper is to examine the structure adjustment of gas giants' 
envelopes after major merger events rather than a statistical study on the coalescence prob- 
ability. We do not consider collisions in which the intruding embryos pass through the 
envelope of gas giants without losing a significant fraction of their heavy elemental masses. 
For merger models, we choose representative rather than a realistic distribution of models to 
illustrate the possibility and character (rather than probability) of two sets of outcomes, i.e. 
the intruding embryos either disintegrate in the gas giants' envelope or impact onto their 
cores. 

For simplicity, we consider mostly head-on collisions. Impacting planetesimals much 
smaller than the giant planet would likely suffer ablative disintegration during their passage 
through the gaseous envelope. On the other hand sufficiently massive impactors may pene- 
trate deep and reach the core. A projectil e's disintegration becomes li kely when it collides 
with an air-mass comparable to its own (JKorycansky fc Zahnld 120051 ). For head-on colli- 
sions, the distance traveled in a gas giant's envelope before disintegration is approximately 
D ~ {pe/pg)Re where pe and pg are the average density of the embryo and gas giant, and 
Re is the radius of the embryo. Reaching the center of a planet {D > Rp) requires a suf- 
ficiently large or massive impactor {R^/Rp > {pg/p^Y). For Saturn-like gas giant planets 
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{pg = 0.5 — 1.5 g cm~^) hit by condensed rocky embryos (pe = 4 — 6 g cm~'^) reaching 
the core becomes feasible for impactors of ~ M^ or larger. We note that the length of the 
path to the core, the air-mass encountered and hence the embryo mass required to reach it 
all depend on the impact orientation. Further, the depth at which an imp actor's mass and 
energy is deposited is sensitive to the impact orientation (JAnic et al.l 120071 ). 



2.7. Mergers of gas giant planets 

In relatively massive disks, the formation of multiple systems may lead to merger events 
due to the following effects: 

Bl) The emergence of the second-generation gas giants can per turb the orbit of their prede- 
cesso rs and induce them to undergo type III runaway migration (JMassetll2008l : IZhang fc Hamilton 
20071 ). During this phase of rapid orbital evolution, planets' orbits may be dyna mically desta- 
bilize d even though they are surrounded by considerable amount of disk gas (jZhou fc Lin 
20071 . also see g22D. 



B2) After the depletion of the disk gas, long-term dynamical instabilities induced by all the 
gas and ice giants on each other (and the embryos) can lead to their eccentricity excitation 



and orbit crossing (JDuncan et al.l 119871 : IZhou et al. 



20071 ) on a time scale 



log{TjPk 



-5 + 2.2A;n 



(2) 



where Pk is the mean orbital period of the planets, ko ~ 2(g — l)ai/{q + l)r^, q = ai+i/ai 
is the ratio of the semi major axis between planets i and i + 1, and the modified Hill's 
radius r' -^ = iMpjaj) -\- Mp('ajj -i)/3M^)^/^a,; takes into account the mass of two interacting 



proto-planets (IZhou et al.l 1200 



YE). 



B3) Eccentricity growth and orbit crossing widen planets' feeding zones but reduce their 
collisional cross sections to its geometric values. Close encounters can eventually lead to 
velocity dispersions comparable to the planets' surface escape speed Ksc- Interior to the 
region where the planets' Keplerian orbital velocity Vkep ^ Ksc, they und ergo repeated close 
encounters until they eventually collide with each other (JLin fc Idalll997l ). 



3. Models of giant impacts onto gas giant planets 

Having established that GIMs are expected to commonly occur among long-period gas 
giants, we use two methods to simulate giant impacts between gas giants and planetesimals 
and mergers between two planets. In the first method, we model collisions into gas giant 
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planets using smoothed particle hydrodynamics (SPH). It is a Lagrangian particle method 
suitable for modeling highly deformed flows and shocks evolving in three dimensions. Because 
it only tracks material particles of interest, modeling flows into empty space can be done 
economically without the need for including vast volumes of all potentially occupied space as 
in a grid type code. This versatile boundary-free multi-dimensional hydrodynamic scheme 
is well suited for the simulations of colliding astrophysical and planetary bodies. 

A few dynamical times after a giant impact, the gross deformation gravitationally settles 
back to a new roughly symmetrical equilibrium shape. Over short timescales and for densities 
near that of the planet, the SPH model provides a good description of the planetary response. 
We examine the long term thermal adjustment of the gaseous envelope using a ID (spherically 
symmetric) Lagrangian radiative hydrodynamic scheme. 



3.1. 3-D Smoothed particle hydrodynamic models 

SPH is w idely used for modeling flow s in a variety of astrophysical and planetary con- 



texts (see e.g. lMonaghanlll992l : lBeng|l99d . for reviews of the me thod). It has been used for 



a variety of problems including modeling stellar col lisions (e.g 



as impacts between condensed solar system objects (JBenz et al. 



1986 : 



Benz fc Hills 1987) as well 



Love fc Ahrens 1996 



Benz fc Asphaug||l999l : lAsphaug et al.lll998l ). The model we use here employs a hi erarchical 
tree to compute self-gravity and is a variant of the code used in iBenz et al.l ( 1l986l ). 



To simulate impacts into gas giant planets with cores, we must model the hydrodynamic 
response o f both condens ed and gaseous matter. We construct r nodel p l anets using the 
Tillotson ( JTillotsonl 119621 ) equation of state (see Appendix II of iMeloshI Il989l ) with the 
material parameters of iron to simulate the structure of the condensed core. The algebraic 
expression of this equation of state is convenient for use in SPH schemes. 

We treat the giant planet's gaseous envelope using a perfect gas equation of state, i.e. 
Pg = {'-f — l)pu where Pg, p and u are the pressure, density and specific internal energy 
respectively and we have used 7 = 5/3. Gas giant planetary models are constructed by 
surrounding an iron core with a gaseous envelope and letting the system relax to a state of 
hydrostatic equilibrium. In regions with modest density, the smoothing length is adjusted 
to maintain a specific number of neighbors. For numerical convenience, we employ a lower 
limit on the density of pi = 0.5gcm~^ such that the expansion of the smoothing length is 
suppressed for particles with density below this value. This density "floor" prevents gas 
in the individual fragments (or the envelope) from dispersing and diffusing out to tenuous 
medium. Since there is limited envelope expansion and very little mass loss associated with 
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models presented here, this numerical prescription does not affect our results. 

Finally, we model large planetary impactors using the Tillotson material parameters 
for both iron and basalt. In principle, impactors might be primarily composed of water 
ice. Impactors' composition for a given mass can determine whether they disintegrate or 
impact onto the core. As we are simply interested in illustrating the gross characteristics 
of merger events (rather than numerically measuring the critical impactor mass or impact 
parameter which determine specific outcomes) our main conclusions are not strongly tied to 
the particular material properties chosen (e.g. basalt, iron, ices) or the use of head-on vs the 
more common oblique collisions for particular impact models. 

As discussed in ^ impactors into gas giant planets range from planetesimals and em- 
bryos to full grown planets. To explore a few of these possibilities, we have performed three 
SPH simulations of impacts into a 'Saturn-like' gas giant planet of mass Mt = lOOM^ which 
includes an iron core with a mass Mt,c = lOM^. The first and second models are designed 
to consider the effects of giant impacts onto a gas giant planet by a planetary impactor with 
Mj = lOM^ and Mj = 25M^ respectively. The impactor is assumed to be entirely made of 
basalt such that its core mass M/ ^ = Mj. A third model examines the merger of two gas 
giants each with Mt = lOOM© and Mt,c = lOM®. 

For long-period planets, it's not likely for the high-velocity collisions to occur. Therefore, 
we set the impact speed fimp = fesc, where the two-body escape velocity is fgsc = {2G{Mi + 
Mt)/{Ri + Rt)Y^'^ for the models in this paper. Model parameters and results of these 
simulations are summarized as "S" series in Table [H 

As expected in relatively low-speed, parabolic collisions, the impactor always merges 
with and is consumed by the gas giants. The gross morphology of the collision product 
and the delivery of the impactor mass and energy into the giant planet are illustrated in 
snapshots of the Sb model shown in Figure [H 

During the initial stages of collision the impactor excavates a crater in the gas giant's 
envelope. Over a few dynamical times the initial deformation propagates throughout the 
planet the planet and body begins oscillating globally and relaxing as the impactor travels to 
planet's center. In the course of this passage, the impactor is flattened and sheds mass before 
eventually reaching the core. As discussed above, a 10-M® impactor is sufficiently massive 
that a significant fraction of its initial mass reaches the core at once. The impactor then 
equilibrates with the core, material settles back into a new state of hydrostatic equilibrium 
and the planet assumes the expected larger spherically symmetric shape. 

A portion of the impactor and impact energy are deposited high up in the gaseous 
envelope. Eventually sedimentation of this mass may release additional energy on a much 
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Model 


Mt 


Mt,c 


Ml 


M,,, 


Mf,c 


M/,, 


Rf/Ri 




(Me) 


(Me) 


(Me) 


(Me) 


(Me) 


(Me) 




Sa 


100 


10 


25 


25 


35 


90 


1.08 


Sb 


100 


10 


10 


10 


20 


90 


1.06 


Sc 


100 


10 


100 


10 


20 


180 


1.25 


La 


100 


10 


25 


25 


35 


90 


1.02 


Lb 


100 


10 


10 


10 


20 


90 


0.996 


Lc 


100 


10 


100 


10 


20 


180 


1.59 



Table 1: Simulations of collisions between the Saturn-like gas giant and embryos with masses 
of 25 Me and 10 Me, and another identical Saturn- mass gas giant, which are denoted as 
model a, b and c respectively. The models named as "S" correspond to the SPH calculations, 
while those named as "L" correspond to the one dimensional LHD calculations. All the 
models are for parabolic, head-on, low-speed collisions. For the composition of the impactors, 
the embryos are made of Tillotson basalt, and the gas giant contains a 10-Me iron core and 
a 90-Me gaseous envelope. Quantities Mt, Mt^c, M/, Mj^c, Mf,ci and Mf^g correspond to the 
target' total mass, its core mass, impactor's total mass, its core mass, the planet's final mass 
in the core and gaseous envelope respectively. Planets' initial and asymptotic radii when a 
quasi hydrostatic equilibrium is re-established after the GIMs are represented by i?j, and Rf 
respectively. 
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Fig. 1. — Collision between a Saturn-like gas giant and a 10-M^ impactor (Run Sb). Gas 
particles are shown as grey. Core particles are overplotted in red, and the basalt impactor 
is overplotted in blue. 
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longer time scale. While the precise partitioning depends on impactor mass, velocity and 
orientation, it is clear that large impacts can deposit mass and energy in large parcels into 
both the core and envelope. 

In addition, the final volumetric radius of the SPH result is in good agreement with 
those of our one-dimensional LHD model (see below) suggesting that longer term thermal 
evolution may be encapsulated using a more computationally convenient approach. 



3.2. The ID Lagrangian radiative hydrodynamic (LHD) impact model 

The 3D-SPH scheme is not ideally suited for the demanding computation of the long- 
term evolution of the gas giant planets. The results generated from the preliminary SPH 
scheme indicate that internal flow quickly isotropizes the density distribution of the merger 
product (see below). In the limit of negligible rotation, the structural response and evolution 
of the envelope, due to the release of thermal energy during the core coalescence, can be 
computed with a ID Lagrangian hydr odynamic (LHD) scheme. This approach was firstly 



adopted by Wuchterl ( iWuchterll Il99ll ) in his calculation of the planet formation problem. 
Although a similar approach is used, our independently developed LHD scheme is based 
on a more conventional prescription for the convective heat transfer and an alternative 
approximation for the computation of the thermal evolution of planets (see §3.2.2p . 



3.2.1. Governing equations for the LHD scheme 

The basic equations for the LHD scheme include three conservation laws of momentum, 
mass and energy, and the equation of state for the relationship among the thermodynamical 
quantities. When writing them in the Lagrangian scheme, these equations are as follows: 



(3) 
(4) 
(5) 



du 

dt 


1 
Po 


(R{r,t)V d{P + q) a^ 
\ r J dr dr 

OR 




1 


1 fR{r,t)ydR 




P 


Po V r J dr 
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" (^)l-^-^ 



dt 



e = f{P,p) (7) 

where the symbols denotes their usual meanings. Note that r denotes the element coordinate 
at some particular time which corresponds to the fluid density po- Usually, these two values 
are viewed as our initial condition. While R{r, t) means the coordinate at time t for the 
same element which initially located at r. In order to smooth the shock, we included the 
artiflcial viscosity denoted as q. Dissipation associated with the artiflcial viscosity is added 
to equation ([6]) to ensure the conservation of total (kinetic, potential, and internal) energy. 
Since we are primarily interested in cohesive collisions and mergers with no mass loss, the 
gravitational potential includes only planet's contribution. 

When discretizing the equations, we use a staggered grid in space. As the code is written 
in Lagrangian scheme, the coordinates will change with the gas flow. Thus, the interval of 
the grids may become very small (for example, in the piling up region) or become very large 
(when the gas is heated by the giant impact and expands). That would not only reduce 
the time step largely, but also make the differentiating between the grid inaccurate. The 
solution is to readjust the mesh after every time step, i.e. merging the nearby two grids 
when the interval becomes small, and splitting one grid into two when the interval becomes 
unacceptably large. 

During the initial evolution after each major impact, we solve the continuity and mo- 
mentum equations explicitly. The time steps for the numerical integration is chosen appro- 
priately so that the Courant condition is satisfled. In order to avoid numerical instability, 
we implicitly solve for the energy equation, which can be converted to the form: 

dT P6 dp 1^ -^ 
at p'^a at p 



3.2.2. Thermal energy transport 

In the SPH scheme, the equation of state for the gas is approximated by a 7 = 5/3 
polytrope and also by an ideal gas equation of state. Kinetic energy is dissipated into 
thermal energy during the impact. Subsequently, the thermal energy is transported to the 
surface and the merger product contracts on a Kelvin-Helmholtz time scale. 
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One advantage of the LHD scheme is its versatihty which allows us to incorporate 
various prescriptions for the equation of state and consider the influence of phase transition, 
elemental separation, density and molecular weight distribution. These effects determine 



wheth er convection may be stabilized in some regions of the envelope. The ISaumon et al. 



( I1995I ) equation of state is used, which especially takes into account the non-ideal effects of 
hydrogen at high pressure. 

For the transported energy flux F in Eq. ([2]), we considered both radiative and convective 
energy transfer. In t he radiative reg i ons, F is calculated using the radiative diffusion equation 
with a flux limiter (JBurkert et al.l |2005| ) . For opacity, we include the contribution from 
grains which is generally sev eral orders of magnitude larger than that from gas molecules 
(jAlexander fc Fergusonlll994l ). 



In our LHD, convection is calculated using the standard mixing length prescription 
(IKippenhahn fc WeigertI 119901) . In the quasi hydrostatic planetary structure models devel- 



oped by 



Pollack et al. 



( 1l996l ). convection is assumed to be efficient and the gas attains an 
adiabatic state throughout the envelope. This assumption is adequate in the envelope of 
present-day gas giants where the density is sufficiently high for convection to be effective. 
However, in the simulation of giant impacts, the envelope may expand greatly to attain a 
low density state in which convection is unlikely to be efficient (see Figure [2]). Therefore, we 
adopt the more realistic mixing-length treatment for the convection during both expansion 
and quasi steady contraction phases. This effect is important in determining the rate of 
radiation transfer, the extent of envelope expansion, and the rate of planetary contraction. 



In contrast to the prescription developed by IWuchterll (1199 ll ). we apply the standard 
mixing length prescription under the assumption that the background gas is in a quasi- 
hydrostatic equilibrium state. This approximation is adequate for the stage during which 
the magnitude of the radial velocity of the gas {\u\) is smaller than the convective speed 
Vcon- However, during the dynamical transition phase after a giant impact, the expansion of 
the envelope is supersonic over some regions of the envelope and the efficiency of convection 
may fall below that estimated with the mixing length prescription. Although the validity 
of the mixing length prescription is weakened, it does provide an upper limit on the rate of 
energy loss and a lower limit on the extent of the envelope's expansion. 



3.2.3. Initial parameters for the LHD models 

We assume the planet is initially in a hydrostatic equilibrium. In order to construct 
an initial model for a gas giant planet, we use Runge-Kutta integration method to generate 
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the internal structure of the planet. In these calculations, constant entropy is assumed to 
represent efficient heat transport in the convection zone. 

The initial model is relevant for a planet with a total mass Mp = 100 M^ and radius 
Rp = 6.2 X 10^ cm. It contains a core with Mcore = 10 M^ and -Rcore = 2.2 R^. The surface 
temperature of this planet is set to be Te = 150 K. Such a model corresponds to a long-period 
gas giant planet with an equilibrium temperature at around 5 AU from a solar-type star. 
These initial conditions are similar to those used in the SPH simulations. 

Starting with these initial conditions, we further evolve the model for a few years and 
compare our results with those obtained by planeta ry structure and evo lution method (un- 



der the assumption of hydrostatic equilibrium, cf. iPoUack et al.l ( 119961 )). This computed 
duration is long compared to the dynamical time scale but short compared with Kelvin 
Helmholtz time scale. Although the results obtained with these two prescription agree well 
in the dense inner envelope, they differ slightly in the outer tenuous regions. The main cause 
for this difference is that convection is unlikely to be efficient in the low density outer regions. 
Although such an approach is more time consuming, we adopt a more realistic mixing length 
prescription for our LHD models. 



4. LHD simulations of Giant Impact and Merger of Gas Giants 

The results of SPH simulations show complex flow pattern during giant impacts of 
embryos with gas giants and during mergers of two gas giant planets. In a few dynamical 
times, the structure of the core and gaseous envelope quickly restores to spherical symmetry 
during head-on collisions. In this limit, we utilize the LHD scheme to study the evolution of 
the GIM products' structure over a larger dynamical range. 

In this series of simulations, we adopt the same model parameters as we used to simulate 
head-on parabolic collisions with the SPH scheme. As expected (and illustrated in §3.ip . 
after initial transitions, spherical symmetry is rapidly re-established. SPH simulations also 
indicate that relatively large proto-planetary cores can survive passage through the gas 
giant's envelope and merge with the core. Based on these results, our 1-D LHD is sufficient 
to simulate the head-on GIMs. 

The results of three sets of simulations are summarized in Table [T] marked as "L". They 
are three sets of models in which the impactor is a 25 M® proto-planetary core, a 10 M® 
proto-planetary core, and a Saturn-mass gas giant. All of these mod els are assumed to be in 



hydro static equilibrium prior to the impact. From virial theorem ( iKippenhahn &: Weigert 



1990) we find the total energy for a polytropic equation of state is 
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W = Eg + E, = {(- 1)EJC (9) 

where Eg and Ei are gravitational binding and internal energy, C, = 8(7 — 1). For 7 = 
5/3, Ei = —Eg/2. At the onset of the simulation, we impose a burst of thermal energy, 
Ea = Mjvf /2, in order to take into account the energy released in the vicinity of the core. 
Therefore, the total energy becomes W = Ea — Eg/2. We also adopt the assumption that 
the cores and embryos are completely merged during the head-on impacts. This assumption 
is supported by the results of SPH models Sa, Sb, and Sc. 

For a impactor with mass of Mj = 25M^, a parabolic collision (model La) leads to 
a merger without any loss of the gaseous envelope. The gas giant's photospheric radius 
expanded slightly. Analogous dynamical parameters are imposed for a 10 M^ impactor. In 
this case, the gaseous envelope is heated by the impact energy and the interior temperature is 
increased by a order of magnitude. However, the gravity of the gas giant is also increased due 
to the accretion of a 10 M® proto-planetary core. Therefore, the expansion of the envelope 
due to the impact heating is suppressed and the gas giant contracted a little bit for this 
model. In Model Lb, the asymptotic radius after the impact is in fact slightly smaller than 
its initial radius as a consequence of the increase in the core mass. In the parabolic model 
Lc, the merger product produces a much larger photospheric radius than the gas giants' 
initial size. Nevertheless, the entire envelope is retained due to a comparable increase in the 
merger's thermal and gravitational binding energy. 



4.1. Mixing of cores and envelope material 

In the discussion above and in §3.1^ we showed with SPH that the impactors with masses 
of 10 Mq and 25 M^ had adequate mass and intruding velocity to reach the core of the target 
planet without total disruption. In this subsection, we consider the mixing between the two 
cores and that between the core and the envelope. The efficiency of mixing determines the 
evolution of the merger's internal structure and mass-radius relation. 

Massive impactors have M/^c's which are significant fractions of or comparable to the 
Mt^c- The dissipated energy Ea is comparable to the gravitational binding energy of the 
target gas giant. If this energy is efficiently converted to thermal energy throughout the 
merged core, temperature in the core would increase by 

ATcore ^ (7 - l){fl/Rg){Ea/Mf,,) ~ {Vi^Jv,spn{M,jMf,,) (10) 

where Tj is the core temperature prior to the impact. 



-20- 



In the models La, ATcore > 5 x 10^ K which is most hkely to be larger than that (Tevap) 
for heavy elements to fully evaporate into a gas phase. In this limit, the vapors of heavy 
elements would thoroughly mix with the nearby hydrogen gas in the envelope. In the SPH 



mode ls where the core and the embryos are approximated as fluids with Tillotson ( iTillotson 



19621 ) equation of state for iron and basalt. 



At lower temperature, the mixture of composite elements undergoes phase separation. 



Altho ugh the phase separation diagram is well established for H-He mixture (jStevenson 



19761 ). its analogs for H-Fe or H-Basalt mixture are not available. Nonetheless, we expect 
the heavy elements in the gas phase to separate from hydrogen into a two-phase medium as 
the temperature decreases below a critical value T2p. During the ph ase separation, t here is 



an exchange between Gibb's free energy and the internal energy Ei fIStevensonI Il976l ) . 



After phase separation, density contrast between the two phase increases with additional 
reduction in temperature. If the heavy elements can assemble liquid drops centered around 
the nucleation sites and grow to sufficiently large sizes, they would precipitate and settle to 
the core (Nucleation is a poorly understood critical phenomenon and it can proceed over 
a considerable time scale). Further reduction in planet's internal tem perature below some 



threshold Tcond would lead to the condensation of the heavy elements (lHubbardlll984l ). 



4.2. Core growth after impactor's disintegration in the envelope 

We consider that for less massive {Mj < a few Mq^) impactors, especially for oblique 
collisions, the impactors may disintegrate before reaching the core. The distribution of heavy 
elemental deposition dMi^ep/dr depends on the impact angle, internal density inside the gas 
giant, and the mass of the impactor. Without the loss of generality, we adopt a simple 
algebraic expression such that the rate of density added at each radius is 

CPdepv / Pnor / ^ / ' ' dep \ \ / ''imp \ /i i \ 

— Ill t- 1 1 gxp( -) (11) 



dt Tdcp V V ^'"dcp // Tdep 

where timp and r^ep represent the epoch of the impact and the duration of disintegration, 
Tdep and Arrep correspond to the main and spread in the location of disintegration, 77 is a 
power index which describes the confinement of disintegration, and pnor is a normalization 
factor for matching the impact's total mass. The actual value of r^ep is determined by 
the disintegration condition, which may have a wide dispersion due to different impact 
trajectories and coUisional speeds even for the embryos with similar masses. 

After the disintegration is completed, dM^ep/dr = in J r'^{dpdcp{r)/dt)dt. The associ- 
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ated energy dissipation rate is 

drdt r drdt 

and change in the local temperature is 

where 

fJ'' = {Hdep{dMdep/dr) + /ienv47rpr2)/(Mdep/c/r + iTipP) (14) 

is the modified molecular weight. 

Depending on the model parameters, the magnitude of AT may not be adequate for 
the new temperature to reach T2p, in which case, the heavy elements are retained their own 
identifies in the form of emulsions. Around the nucleation site, heavy elements congregate 
and form droplets. Due to the buoyancy effect, these droplets (with density p^ greater than 
the local gas density Pg) have a tendency to precipitate. In an convective region where the 
local convection speed is fconv, droplets with sizes s larger than a critical size 

Ssettlc ~ PgVconv/pd9 (15) 

(where g is the local gravity), would sediment with a speed 

^settle ~ (s/Sscttlc)^^^^^conv ~ {sQPd/ Pgf^"^ ■ (16) 

Smaller droplets may also be dragged by the convective eddies, diffuse to and become part 
of the core. 

The actual size distribution of the droplets may depend on not only how the impactors 
break up, but also how their fragments subsequently interact with the ambient gas in the 
envelope. Equation flTTl) is based on the assumption that the impactor is completely disinte- 
grated down to the molecular scales which then form liquid droplets around the nucleation 
sites. These liquid droplets coagulate as they collide through dispersive motion or differen- 
tial settling. They also fragment due to both collisional fragmentation as well as Kelvin- 
Helmholtz and Rayleigh- Taylor instabilities at their interfaces with the envelope gas. For 
sufficiently small droplets, these instabilities are stabilized by their surface tension. Subli- 
mation (from liquid to gas phase), condensation (from liquid to condensed state), chemical 
reaction, and sedimentation can also contribute to the size distribution. 

Despite all these uncertainties, we consider a general power-law size distribution in 
which dN/ds oc s'^. For solid particles in a gas free environment, collisional equilibrium 



is est ablish with A ~_— 3 -5 (IDohnanvi 



19961 : iTanaka et al.lll996l : iKenyon 



2002 
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969l: IWilliams fc Wetherilll ll994J : iFarinella fc Davis 



Weidenschillingl 12004 ) . 



The SPH models presented in §3.11 show differentiation occurs rapidly. This evolution 
is possible if the heavy elements can concentrate into large coherent clouds. However, these 
results may be an artefact of the SPH scheme which uses a separate set of independent 
"fluid particles" to represent the hydrogen gas. This scheme does not take into account the 
Gibbs free energy associated with the mixing of two species. On the scale of the individual 
fluid elements, they are not subject to additional fragmentation. Since this scale is >> Ssettie, 
sedimentation speed becomes a significant fraction of f esc and the time scale for differentiation 
becomes comparable to that for dynamical collapse. 

In principle, the sedimentation of the disintegrated debris should be analyzed with a 
two fluid approach. However, if the fragments are sufficiently small to be well coupled to the 
gas, the flow can treat the mixture as one gas with a modified molecular weight //'. With 
equations ( ITTj) . ( !T2|) . and (JH]), we include the effect of energy deposition, mass loading, and 
molecular weight increases in the LHD models. 

In order to simulate an event in which an embryo disintegrates in the planet's envelope, 
we introduce a new model Lbl. All the initial and boundary conditions in this model are 
identical to those of model Lb, except we set r^ep = '^rcore, Ardep = O.Srdep, and Tdep = lO^s 
in equation f lTT]) . The normalization factor Pnormai is adjusted so that the total mass of the 
impactor remains to be 10 M^. The associated Ea = 2Ax 10^^ ergs which is similar to that 
(2.22 X 10^^ ergs) for model Lb. 

In model Lbl, debris of the impactor is deposited well outside the target gas giant's 
core. Over a few dynamical time scales, a quasi hydrostatic equilibrium is established in 
the envelope and the flow speed everywhere inside the photosphere decreases below 0.1 
Cg. Thereafter planet's interior structure evolve through thermal adjustments as heat is 
transferred from its interior and lost near its photosphere. The efficiency of heat transfer is 
determined by not only radiative diffusion but also thermal convection. 

According to the standard Schwarzschild criterion, the chemically homogeneous plane- 
tary envelope is generally unstable against convection. However, heavy element and energy 
deposition after a major impact modify the local thermal and compositional distribution. 
According to the Ledoux criterion, the flows can be stabilized by the molecular weight gra- 
dient if 



Vt < Vad + (<^/5)V, 



(17) 



where the radiative Vt = {P/T){dT /dP) = {?,KF/AacT^){P/g), the adiabatic V ad 
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(7 — l)/7 ~ 0.4, chemical V^ = {d\nfi/dlnP)s. The subscript s refers to the surround- 
i ng medium. For an ideal gas , both ip{= dlnp/dlnfi) and 6{= —dlnp/dlnT) equal to unity 



( iKippenhahn &: Weigertlll990[ ). 



The initial density distribution for model Lbl is plotted on the left panel of Figur^ 
Along with the density distribution, overplotted are the Vt, and V^ distribution shortly 
after a quasi hydrostatic equilibrium is reached. At this epoch, there are three interesting 
regions. 

a) Outside r = 3.3 x 10^ cm, the envelope is convectively unstable because both 
Schwarzschild and Ledoux criterion for stability are not satisfied. The impulsive heat depo- 
sition drives vigorous convection to transport heat and to mix the heavy elements with the 
envelope gas. 

b) In the region r = 3 x 10^ — 3.3 x 10^ cm, where the Ledoux criterion for stability is sat- 
isfied but Schwarzschild c riterion is not, the envelope is dynamically stable but vibrationally 



unstable (or over-stable) (IKippenhahn fc WeigertI 119901 ). In this limit, double (thermal and 



chemical) diffusion leads to layered structure in which diffusive radiative layers are sand- 
wiched between convective layers which are well mixed and have adiabatic structures. When 
averaged over many layers, stratification would probably adjust toward a state of marginal 
stability. With an equation of state, dp/p = a{dP/P) — 6{dT/T) + ip^dp/p), the Ledoux 
criterion is marginally satisfied with 

Vtot = Vrad - Vadiabatic " V ^,{i> / 6) ~ 0. (18) 

c) Interior to r = 3 x 10^ cm, the flow is stable according to both criterion and the 
initial convective envelope is stabilized by the heat released and heavy element deposition. 

In both regions b) and c), convection is suppressed with vanishing Vconv After the 
impact, the temperature in these regions is well above Tcond (see Figure [3] panel a) and most 
likely above T2p. Despite the stabilization of regions b) and c) against convection, heat 
is still transferred to the planet's surface through radiative diffusion. On a much longer 
Kelvin-Helmholtz cooling time scale (~ 10 — 100 Myr), phase separation would lead to the 
formation of liquid droplets when the envelope's temperature declines below T2p, provided 
that the condition allows the conversion of Gibb's free energy to internal energy. 

Droplets grow through condensation and coagulation. While the convective stability is 
maintained, relatively large droplets would sediment towards the core with a terminal speed 
Vsettie ~ {s9Pd/ PgY^"^ ~ (s / RpY^'^Vesc- In this Stoke's range, the sedimentation time scale 
(Tsettie ~ (Rp/sY^'^Tdyn, wlicrc Tdyn ~ {Rp/GMpY'''^ IS the dynamical time scale) is shorter 
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Fig. 2. — Gradients for temperature and molecular weight distributions for models Lbl 
(on the left panel) and Lb2 (on the right panel) after the models have reached quasi-static 
equilibrium. The black solid line shows the temperature gradient. The red dotted line 
corresponds to the molecular weight gradient. The density profiles for both models are 
overplotted on the same plots. The green dashed line represents the density profile of the 
initial model. The blue dash-dotted line shows the density profile at the same epoch as when 
those gradients are plotted. 
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Fig. 3. — P-T diagrams for model Lbl (on the left panel) and Lb2 (on the right panel). 
Overplotted is the melting curve of iron. 
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than the age of mature solar- type star, r*. 

Drag on the small particles (with s < Smfp, where the mean free path Smfp — 5 x 
10^^ cm /(p/lgcm^^) for the ambient hydrogen gas) is in the Epstein range. In this limit, 
Vsettie ~ {sg/cs){pd/Pg) ~ {s/Rp)vesc and Tsettie ~ {Rp/s)Tdyn which may be substantially 
longer than r*. 

The above consideration indicates that the debris deposit of moderate-mass embryos 
which disintegrated in the planetary envelopes not only introduces a metal-enriched layer well 
outside the core but also locally suppress convection. These two effects promotes fragments 
sedimentation and the growth of their cores. 



4.3. Core erosion after impactor's disruption near the core 

It has long been recognized that compositional stabilization is an important process in 
preserving the heavy-el ement-rich cores ag ainst thermal convection of the gaseous envelopes 



in isolated gas giants (iGuillot et al.ll2004l ). If turbulence in the convective envelopes can 
reach the cores, mixing and erosion would reduce Mt^c- However, this diffusion leads to a 
strong negative molecular weight gradient across the core-envelope interface. This gradient 
would quench thermal convective instabilities, suppress the compositional diffusion rate, and 
preserve the core-envelope separation. 

In §3.11 and §11 we showed that in nearly head-on GIMs, sufficiently massive impactors 
can penetrate deeply into the target gas giants' envelopes and reach their cores. For these 
GIMs we show in this subsection, that the impulsive dissipation of impact energy would 
strongly drive thermal convection despite the compositional gradient introduced by the heavy 
element injection and diffusion. 

In order to demonstrate this possibility, we construct model Lb2, in which, debris of 
the impactor is deposited in the proximity of the target gas giant's core. All the initial and 
boundary conditions in this model is identical to those of model Lb, except we set r^ep = ^core, 
Ardep = 0.25rdep, and r^ep = 1 x 10^ s in equation (fTTj) . The normalization factor pnormai 
is adjusted so that the total mass of the impactor remains to be 10 M®. The associated 
Ea = lMx 10^1 ergs. 

Similar to model Lbl, a quasi hydrostatic equilibrium is established in the envelope and 
the flow speed everywhere decreases below 0.1 Cs after a few dynamical time scale. We plot 
Pg, Vr and V^ distribution at this epoch for model Lb2 on the right panel of Figure [2j By 
this stage, a fraction of the debris of the disintegrated impactors has reached r = 4 x lO^cm. 
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Within this radius, both Schwarzschild and Ledoux criterion for stabihty are violated and the 
heat release near the core strongly enhances the tendency for thermal convection throughout 
the envelope, including the boundary layer region. 

The subsequent evolution of the heavy elements depends on whether the collisions lead 
to a total phase transition of the core. Let us first consider the possibility of a hot core 
with an initial temperature Tcore > ^cond so that the heavy elements are in a liquid or gas 
state. Provided the core-envelope interface is stabilized by an initial compositional gradient 
in accordance with the Ledoux criterion, this differentiated structure would be preserved. 

After the impact, vigorous and efficient convection would induce the interior of the 
planet to evolve towards a new constant-entropy state. During that evolution, if a fraction 
of the dissipated energy fEa is transported into the central region, the fractional change in 
the core would be ATcore/T'core ~ fEa/Ei (In the limit Ea > Ei, ATcore would be given by 
equation [To!) . The time scale for thermal adjustment towards this state is approximately 

Tad ~ fEjinRl^^F,,,, (19) 

where -Fcore is the rate of heat transport into the core (In the limit of high heat conduction 
rate, Fcore would be the convection flux near the core radius -Rcore and it can be computed 
using standard mixing length prescription). During this time, the amount of heat loss from 
the planetary surface is ~ inRpFg-^^ where F^ur is the planet's intrinsic heat flux at its surface 
(Note Fsur does not include the reprocessed stellar irradiation). From energy conservation, 

we find / ~ (1 + {F,^r/F,on){Rp/ Rcore?- 

With sufficient input energy Ea, the core temperature can increase above Tgvap so that 
the heavy elements are evaporated. In the limit of efficient convection, the envelope's tem- 
perature distribution adjusts to establish constant entropy state. Heavy element would mix 
with the gas in regions where the local temperature is above T2p. If this radius extends 
well beyond the initial Rcore, a substantial fraction of the core would be eroded and mixed 
with the envelope. However, if after the impact, T^ore = ^corc + ATcore < T^vap or the region 
with T > T2p is confined to small radius, the core structure would be mostly preserved even 
though convection may be able to penetrate through the initial core-envelope interface. 

We now consider the possibility of GIMs onto gas giants with cool (Tcore < ^cond) cores. 
With sufficiently massive and energetic impactor, head-on collisions lead to T^orc > ^evap- The 
subsequent evolution is similar to that discussed above, albeit a modest fraction of internal 
energy must be used to compensate for the latent heat during the melting and evaporation 
of the core. Less energetic impacts may not be able to evaporate the core, in which case, 
the temperature distribution of the core would be determined by the energy equation ([8]) 
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including heat transport by conduction. Provided the new T > Tcond at the core-envelope 
interface, the mohen heavy elements would be entrained into the convective eddies on the 
growth time scale for Kelvin-Helmholtz instability. This rate may be limited by the poorly 
determined material properties such as surface tension and phase separation. However, if 
T > T2p at the core-envelope interface, mixing between the heavy elements and hydrogen 
may be much more efficient. Since convection in the envelope extents to the interface, the 
mixed gas is redistributed throughout the convective region with T > T2p. 



5. Summary and Discussions 

The main objective of this paper is to suggest that giant impacts and merger process may 
have led to different core-envelope structure between Jupiter and Saturn and the super-solar 
metallicity in their gaseous envelopes. 

We have highlighted several avenues of GIMs by embryos onto gas giants. The sup- 
ply of these large solid building blocks include: in situ embryos during the final phase of 
gas giant formation, and terrestrial bodies cleared by sweeping secular resonances and dy- 
namical instabilities. There are also several processes which can lead to the merger of gas 
giants, including orbit crossing triggered by run-away migration, and long-term dynamical 
instability. 

At their present-day orbits, the velocity dispersion of any residual embryos or plan- 
etesimals is less than the gas giants' surface escape speed. The GIMs between any residual 
planetesimals, embryos, or hypothetic gas giants and long-period gas giants are likely to be 
parabolic and attain limited amount of energy. This modest velocity dispersion also enlarges 
gas giants' gravitational cross section. 

The discussions above indicate that embryos with a mass of several M^ can penetrate 
through the envelope and reach the core of gas giants through head-on collisions. If these 
collisions are sufficiently energetic, they would lead to significant core erosion. Heavy ele- 
ments would be well mixed with the hydrogen in the envelope down to a molecular level. The 
average density and metallicity of the gas in the envelope would be enhanced. Although a 
relatively large amount of kinetic energy is dissipated into internal energy, the heat released 
from this type of GIM is efficiently transported to the planet's surface and radiated away 
because the envelope is fully convective. 

Low-mass impactors (M/ ^ M^ are likely to disintegrate in gas giants' envelope. Energy 
deposition in the planet's envelope suppresses convection interior to r^ep- Heavy element 
deposition also increases the opacity. After the planet has adjusted to a new hydrostatic 
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equilibrium, the impact dissipation energy is more effectively retained by the reduced heat 
ffux through the envelope. Suppression of convection also leads to sedimentation of modest 
to large size embryos. In the limit that a collisional equilibrium can be established with 
most of the debris' mass being in the forms of large particles, the core attains, from the 
post-impact differentiation, a mass su bstantially larger than the critical m ass (a few M^) 



for the onset of efficient gas accretion (jPoUack et al.lll996t llda fc Linll2004j ). 



We suggest these two limiting cases introduce diversity in the core-envelope structure. 
They may also lead to different cooling rates during the subsequent thermal evolution. For 
the core-impact case, the dredge-up of heavy elements does not suppress convection. During 
the cooling process, phase separation and phase transition lead to changes in the Gibbs free 
energy which may lead to the release of additional sources of internal energy to reduce the 
planet's contraction rate. For the envelope-disintegrate c ase, small grains or droplets are the 



dominant sources of opacity (Ide Pater fc Lissauen I2OOII ) and their suspension or depletion 



determines the contraction rate of the gas giant. Analysis for the long term thermal evolution 
of planets will be considered and presented elsewhere. 

We suggest that Jupiter may have been struck by an energetic head-on GIM which 
penetrated deeply into its envelope and lead to its core erosion. With a larger (> lOM®) 
initial core, it would also partly resolve a long standing issue on how Jupiter might acquire 
its massive gaseous envelope with such a small core on the multi-Myr time scale over which 
their natal solar nebula is depleted. 

In contrast, Saturn may have acquired its present big core through GIMs with less ener- 
getic and smaller residual embryos. A temporary stabilization of Saturn's envelope against 
convection may have also prolonged its contraction time scale and enabled the precipitation 
of helium dropout in its envelope. We note that both processes lead to the enrichment of 
metallicity in the envelopes of Jupiter and Saturn. 

GIMs are most likely to occur within the first few Myr of their life span. Energy 
dissipation from such events is likely to increase the intrinsic luminosity and radius of these 
planets. Direct imaging and rare transit observations may provide observational verifications 
that the internal structure of typical gas giants may have been modified by GIMs. We shall 
discuss the GIM's influence on the mass-radius-luminosity of the close-in planets in Paper 
II. 

We have shown possible mechanisms for the dredge up of heavy elements from the 
core and for the promotion of core growth. We have not considered here the long-term 
sustainability of the stir-up elemental distribution. The follow up evolution of the impacted 
gas giant's radius depends sensitively on their cooling efficiency and elemental segregation. 
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Althou gh there are models for contraction and diff erentiation of initially homogeneous gas 



giants (JBodenheimer et al.ll2003l : iBaraffe et al.ll2008l ). the long-term post GIM evolution still 



needs to be investigated. 
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